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We reconsider the dispersive evaluation of the weak matrix elements 
((7T7r)/ = 2| 27,81-^°) in the chiral limit. The perturbative matching is accom- 
plished fully within the scheme dependence used in the two loop weak OPE 
calculations. The effects of dimension eight (and higher dimension) opera- 
tors are fully accounted for. We perform a numerical determination of the 
weak matrix elements using our dispersive sum rules fortified by constraints 
from the classical chiral sum rules. A careful assessment of the attendant 
uncertainties is given. 



Typeset using REVTeX 



I. INTRODUCTION 



In this paper, we present an updated and improved procedure for obtaining analytical 
expressions and numerical evaluations of the matrix elements ((7T7r)/ =2 | Q-j^\K°) in the chiral 
limit. Recall |l| that the ratio e'/e can be expressed numerically in terms of operator matrix 
elements as evaluated at the scale /x = 2 GeV in the MS'-NDR renormalization scheme [0] 
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where 



(Q 6 )(°) = ((7T7r) /=0 |Q 6 |K°) and <Q 8 >< 2 > = <(7r7r) 7=2 | Q 8 \K°) 
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and all other notation is as in Ref. 0. In a previous work ||, it was shown that in the chiral 
limit the matrix element (Qs)n (and also {Qj)y) is expressible in terms of certain vacuum 
matrix elements^ 
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where F^ is the pion decay constant evaluated in the chiral limit and the operators Q\ t , 
are defined as 
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In the above, q = u,d, s, r 3 is a Pauli (flavor) matrix, {A a } are the Gell Mann color matrices 
and the subscripts on Ox, Os refer to the color carried by their currents. Donoghue and 
Golowich showed in Ref. how to obtain dispersive sum rules (which we shall refer to as 
DG1 and DG2 in this paper) for the vacuum matrix elements of the dimension six operators 
Ox and Os- In a later work Q], the presence of higher dimension operators (i.e. those having 
dimension d > 6) was identified and their impact discussed. 

This paper will extend previous work in several significant respects: 



1 In this paper, we work with the operator Qg = s a T^db yubT^u a — \dbT^d a — \sbT^s a j to be con- 
trasted with Qg 3//2 ^ = s a Y^db (ubT^u a — dbT^da^j + s a Ti j UbUbT^d a as used in Ref. ||. In particular, 
one has {(tttt) I=2 \Q 8 3/2) \K°) = 2<(vr7r) 7=2 | Q 8 |^°> • Throughout we define Tt = 7^(1 + 75). 
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FIG. 1. ALEPH determination of Ap(s) 



1. We provide a two-loop determination of the dimension-six contributions to the operator 
product expansion (OPE) for the isospin correlator AIT(g 2 ), 

i J d 4 x e^{0\T(Vr(x)V^(0)-A^x)A^O)) |0> 

= W - gV )An(g 2 ) - q^%{q 2 ) • (5) 

We give our results in MS renormalization, using both NDR and HV schemes for 75, 
and adopting the same evanescent operators scheme as in Refs. 

2. We extend previous work on dispersive sum rules 0] with an updated derivation which 
takes into account the NLO renormalization scheme dependence, and incorporates 
contributions from so-called higher- dimensional operators @. 

3. We perform a dispersive evaluation of (Oi)^ and (P^)^. Our numerical analysis unifies 
input from experiment (the existing data base for the spectral function Ap(s)) with 
rigorous theoretical constraints embodied by the Weinberg and pion mass differ- 
ence H chiral sum rules. An important consequence is the assignment of realistic 
uncertainties to our results. 



II. THEORETICAL ANALYSIS 

In Ref. it was shown how to relate the vacuum matrix elements of 0\ and O s to the 
correlator AIT(g 2 ). In this section we describe an improved version of the theoretical analysis 
of Ref. ||, beginning with some considerations on the central object of the analysis, the 
vacuum correlator AU(q 2 ), defined in Eq. (|5|). Our work exploits a combination of the two 
nonperturbative representations available for the correlator, the dispersive representation 
and the OPE representation at large spacelike momenta. 

The dispersive representation reads 
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where Q 2 = — q 2 is the variable for spacelike momenta and the difference of spectral functions 
is denoted by Ap(s) = [pv,z ~ Pa,z\ ( s )- The ALEPH || determination of Ap(s) for 4m 2 — 
s < m 2 is displayed in Fig. [I]. 

For Q 2 3> Aq CD , AU(Q 2 ) can be represented via the OPE. Through 0(a 2 ) one has 
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(d = 2,4,6,8,10,...) 
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where ad(p) and bd(p) are combinations of vacuum expectation values of local operators of 
dimension d. We list some properties of the OPE: 

1. For d < 6, the coefficients a^, bd are known to 0(a 2 ) for d = 2 and 0(a s ) for d = 4, 6. 

2. For d < 6, a^, 6^ are 0(m q ) or 0(m 2 ) and thus vanish in the chiral limit. 

3. a6 and &6 are related to vacuum matrix elements of the operators 0\ and 0%. 



4. For d > 6, ad, bd are partially known (an analysis of d = 8 can be found in Ref. |IJ for 
the vector correlator and Ref. [jlll for the V-A correlator). In this work their detailed 
form is not needed. Hereafter we denote the collective d > 6 contributions to the OPE 
as ATT(Q 2 ), that is 



AU(Q 2 ) ~ £ 



d>6 
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By virtue of item 3, the d = 6 OPE coefficients b G are of special interest. Here, we 
consider the MS renormalization scheme with NDR and HV prescriptions for 75 and the 
evanescent operator basis used in Refs. |5||J. Including 0(a 2 ) terms we find for N c = 3 and 

n f = 3,4, 

a 6 (/i) = 2ir(a.O B ) lt + A 8 (a 2 s 8 )^ + A^O^^ , 

be(ji) =B s { a 2 s O s )^ + B l {a 2 s O l ) ll . (9) 

The coefficients A\, A§ and B\, B$ are displayed in Table H, in terms of their dependence on 
the renormalization scheme and the active number of flavors. We shall present a derivation 
of the above results in Sect. H A| . An additional piece of information which will be needed 



in our analysis is the form of a§(ii) at 0(a s ) in d = 4 — e dimensions, 

o fl (//,e) = 2ir(a s O s )Jl + k s e^ . (10) 
The scheme- dependent coefficient k s can be found in Table |. 
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TABLE I. Collection of coefficients needed in the calculation at various stages. 





Three Active 


Flavors (nj = 3) 
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Scheme 
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HV 


NDR 


HV 


Ai 
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-10/3 


2 


-10/3 


^8 


25/4 


21/4 


205/36 


169/36 


Bi 


8/3 


8/3 


8/3 


8/3 


Bs 


-1 


-1 


-2/3 


-2/3 


Cs 


-1/6 


11/6 


-1/6 


11/6 


k s 


-1/4 


3/4 


-1/4 


3/4 



A. Determination of a§{ij) and b^(fi) 

In this section we derive the expression for the dimension six terms in the OPE for 
AU(Q 2 ) at NLO in QCD. We begin by defining an amplitude M. which describes the coupling 
of isospin vector and axialvector currents via W-boson exchange, 

M^^Jd'x V(x, M 2 W ) (0|T (V 3 M (*) ^(0) - A$(x)A^)) |0) . (11) 

Our strategy is to perform two different analyses of this amplitude at NLO in QCD (one of 
which involves the vacuum correlator we want to study), and then match the two: consis- 
tency then determines a 6 (fj,) and b & (n). 

Within the first approach, we write the amplitude M. in the language of an effective 
theory, involving the operators Oi, Og and their Wilson coefficients ci, c 8 , 

M^^^lcM^. + csi^iOs),} , (12) 

where /i is the renormalization scale. We define the effective theory via dimensional regular- 
ization within MS renormalization. The specification of how 75 is treated in ^-dimensions 
(NDR or HV scheme) and of the evanescent operator basis is needed to uniquely define the 
effective theory. We follow here the prescriptions given in Ref. ||. The Wilson coefficients 
C\{p) and Cg(yu) are found by performing perturbative matching of the full and effective 
theories at scale fi = Mw and then evolving them according to the renormalization group 
equations 

A*T" C M = (lik ~ 2^fj5 tk ) c e (ji) ■ (13) 
dji 

Here {7^} is the anomalous dimension matrix for the operators 0\$, while 7j is the weak 
current anomalous dimension^. For an NLO analysis we need ci$(Mw) up to non-logarithmic 



2 7j is nonzero at two loops in the HV scheme, when subtracting minimally. In this work we stick 
to this version of MS — HV renormalization scheme. 
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terms of order a s , and also the two loop anomalous dimension. The matching calculation 
required here has been performed in Ref. || and results in 



Cl (//~il%) = l + 0{a 
where d s is given by 



si 1 



d. 



tjr \ 3a s (M w ) /3 „ , 
c 8 (/i ~ M w ) = ^— ^ - + 2d, 

07T \Z 



-5/6 (NDR) 



1/6 (HV) 



(14) 
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The anomalous dimension matrix {7^} is parameterized in terms of the number of active 
quark flavors rif and the number of quark colors N c . At next-to- leading order (NLO), it has 
the form 



a s{v) (0) , ( a s(n) \ (i) 

7nlo = — ; — 7 + — ; — 7 



47T 



47T 



(16) 



The next key observation is that the anomalous dimension matrix for the operators Oi$ 
can be inferred by the restriction of the full 10x10 matrix of Ref. |||6| to the Q-j^ subspace. 
The two sets of operators 0\$ and Qi^ have the same Dirac structure, but differ in the 
flavor structure and in the basis for the color structure. Since only current-current diagrams 
contribute to the anomalous dimension matrix of C?i,8j an d these contributions are insensitive 
to flavor, we only need to handle the different colour properties. The difference in the color- 
structure basis is taken care of by a simple linear transformation. Denoting by the 
anomalous dimension matrices of Refs. [|,[| restricted to Qr,s, those we need for Oi )8 are 
given by 



with 
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(17) 



(18) 



Using the above ingredients and the NLO evolution operator we have calculated the 
Wilson coefficients in MS renormalization for both NDR and HV schemes, taking n/ = 3,4 
and N c = 3. Upon using the expression for the 2- loop running a s , we have then expanded 
c i,8(a0 m powers of a s (fi), finding 



Ci(/i) = 1 + 



7T 



3A 1 , M, 
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+ ..., 
(19) 



where A\, A 8 are the scheme-dependent coefficients of Table 
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The alternate analysis of the amplitude M relies on its expression in terms of the corre- 
lator An(Q 2 ): 

M = " G f 2w r dQ 2 ^ Q * An(Q 2 ) • (20) 

32^2^2 Jo Q 2 + Myy y 1 

We partition the amplitude M as 

M = MM + MM (21) 
where the component «M<(/V) arises from contributions with Q < fi, 

MM = 32 ^ 2p2 ff dQ 2 Q 4 AH(Q 2 ) + 0(^/M 2 w ) (22) 

and the component M > (/i) contains the contributions with Q > /i, 

= zGfM^ /■« 2 Q 4 An( Q2 ) (23) 

>VW 32^27r 2 F2V Q2 + M 2 , w; v ; 

Employing the OPE for Q 6 An((5 2 ) in 7Vl>(/i) and performing the Q 2 integration, one can 
compare this representation with Eqs. (|T2])-( |T9|) . Requiring the consistency of the two ap- 
proaches yields the expressions of Eq. ([]) for a 6 (/x) and b 6 ([i), as well as the sum rule DG1 
(cf. Eq. (|25D below), to be derived in a different way in the next section. 



B. The Two DG Sum Rules 

In this section we present a short derivation of the two DG sum rules ||, including 
effects which were neglected in the original derivation. Let us start with the first sum rule. 
Consider the V-A correlator in coordinate space for x — ► 0, regularized dimensionally 

(O l ), = (0|T (^(0)^,3(0) - ^(0)^,3(0)) |0>„ 
(d- l)/i 4 " d 



/ dQ 2 Q d AU(Q 2 ) . (24) 
Jo 



(47r) d / 2 r(d/2) Jo 

Let us now break up the above integral at a scale yU 2 , within the region of applicability of 
perturbative QCD. The integration over Q 2 < fi 2 is UV finite and can be performed in d — 4. 
The integration for Q 2 > fj, 2 generates UV singularities. Using the OPE representation for 
AU(Q 2 ) in this region, one sees that the UV divergence is related to the Q~ 6 term in the 
expansion, while An(Q 2 ) leads to an UV finite term. Employing Eq. (0) for a 6 (/i, e), 
performing the Q 2 integration and subtracting the 1/e pole according to the MS scheme, we 
obtain 

(d), - g^<«.0B> M = U») ( DG1 ) (25) 
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where the scheme- dependent coefficient C 8 can be found in Table HQ and Ji (/i) is defined in 
terms of the dispersive integrals h{^), -Hi(aO as follows 



where 



A (a*) 



Mm) 



(is s In 



(4tt) : 



S + /i 



dQ 2 Q 4 An(Q 2 



rfQ 2 q 4 An(g 2 



Evaluating the OPE of Eq. (|7|) at the point Q = \l yields the second DG sum rule, 
27r(a s G 8 ) M + A 1 {a 2 s 1 ) tl + Asia^Os)^ = 2ica.(ji)I a (n) (DG2) 



where we have made use of Eq. 
where 



and define I 8 (/i) as 
' r / 8 (/i) - H 8 (n) 



ds s z 



10 s + /i 2 

hm = /i 6 An(/i) . 



Ap(s) = // 6 AII(//) 



(26) 

(27) 
(28) 

i 

(29) 

(30) 

(31) 
(32) 



In Ref. the terms ifi i8 (/i), subleading at high /x, were neglected. They encode the effect 
of higher dimensional operators and lead to potentially large effects |fHJlT 
We can summarize the work thus far via the linear relations 



-3C 8 a s (/i)/87r 



\A l a s {^)/2Ti (1 + A 8 a s (fi)/2n) 



(33) 



which allow us to keep track of the scheme dependence of the matrix elements at NLO. As a 
check on our calculation of the coefficients Ai )8 , C 8 , we can derive the relation between the 
HV and NDR matrix elements of 0,7$. To do so, we need to solve Eq. ( |33"D for (Oi^)^ in 
both schemes and then convert to using the matrix M of Eq. (|i~8"l). We find 



fiQrh 



HV 



(Qr) t 



NDR 



in agreement with the result given in Refs. RIT3 . 




NDR 



(34) 



3 Although we have employed three distinct quantities k s ,d s ,Cs to indicate how scheme- 
dependence enters in different parts of the analysis, they are in fact related by C 8 = 1/3 + 2k s = 
3/2 + 24. 
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C. Additional Comments 



Let us take note of the differences of the preceding analysis from that in the previous work 
by two of the present authors 0. In the earlier work, the information about perturbative 



corrections in the NDR scheme was taken from the paper of Lanin et al. JTj] in the literature 
on QCD sum rules. This is now seen to differ in the renormalization conventions and 
use of evanescent operators from the choices made in the two loop analysis of the weak 
Hamiltonian. In the present work we have adopted the conventions of the weak interaction 
studies. One fortunate benefit from this change is that the results are less sensitive to 
higher order perturbative corrections than in the previous analysis. In the earlier work 
we also did not present all results for the HV scheme, whereas this is accomplished in the 
present study. Finally, while the previous work led directly to the uncovering of the effects of 
higher dimension operators, these effects were not explicitly incorporated into the formulas 
or numerical work of that paper. In the present paper we have explicitly included these 
effects. 

The integrals /sOO an d AO 4 ) wm be central to our analysis. They are related by 



J 8 (/i)=/i 2 ^/ 1 ( /U ) , (35) 



and with the aid of the Weinberg sum rules |7]], it is possible to express /sO 4 ) m several 
equivalent forms, e.g. 

I 8 (//) = -// 4 / ds—- 2 Ap(s) . (36) 
Jo s + n z 

Moreover, there is an alternative pathway to Eqs. ([27|) , ([H]) . One simply evaluates two 
contour integrals in the complex s = q 2 plane, 

i ds K k (s, n) An(s) = (A; = 1,8) , (37) 
where the weights Kk{s,fi) are 

-Ki(s,/x) = s 2 ln , K 8 (s,fjL) = — ■ — ■= . (38) 

s s + p z 

The contours r§, I\ are shown in Figs. [|(a),(b), with the understanding that the radius of 
each outer circle is to be taken to infinity. The unitarity cut in Figs. 0(a), (b) is denoted by 
a wiggly line along the positive real axis. The pole in K&(s,fi) at s — — /i 2 is denoted in 
Fig. ^(a) by the cross on the negative s axis. The cut in Ki(s, fi) occurring for — /i 2 < s < 
is denoted by the bold straight line. 

Obtaining numerical values (along with meaningful error estimates) for IsOO and ii(/x) 
turns out to be a highly nontrivial task. It is this problem that we turn to in the next 
section. 
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(a) (b) 
FIG. 2. Integration contours for (a) ig, (b) I\ 



III. NUMERICAL ANALYSIS 

The goal of this section is to arrive at numerical estimates for the input vector on the right 
hand side of Eq. fl3"3"|). We shall accomplish this by means of a constrained evaluation of the 
dispersive integrals I\^{p) using an approach which we call the residual weight approxima- 
tion. This determination can be applied over a range of possible scales, e.g. 2 < p(GeV) < 4. 
In fact, we choose a scale /i = 4 GeV sufficiently large for higher-dimension effects to be 
negligible and then employ a two-loop renormalization group analysis to evolve our result 
down to scale p = 2 GeV. The decision to set #1^(4) = is supported by a preliminary 
study using Finite Energy Sum Rules (FESR) flSfl . Work on this topic continues and will 
be reported on in a separate publication fi~6|| . 

A. Residual Weight Approximation 

Direct evaluation of h(p) and Is(aO y i a Eqs. (^]) and (|3l"l) is impossible since data for the 
spectral function Ap(s) exists only in the interval 4m 2 < s < m 2 . Even direct integration of 
the contribution from the region s < m 2 is potentially problematic since strong cancellations 
significantly enhance the impact of experimental errors. This is compounded by the fact 
that the weights iT li8 (s, p) are strongly increasing with s and hence weight the high-s region, 
where data errors are large, more strongly than the low-s, low error region. As an example, 
defining /i 8 (so,/i) to be versions of Eqs. ( p?D and ( pT| ) containing spectral contributions 
only up to s = sq, the errors on Ii.s(so, /i), determined using the ALEPH covariance matrix, 
already exceed 50% of our final central values for s ~ 2.1 GeV 2 . The possibility of obtaining 
reliable determinations of /i,g(^) thus depends critically on our ability to impose additional 
constraints on these integrals. 

Although Ap(s) is not known above s = m 2 ., and is determined with insufficient precision 
for our purposes between 2 GeV 2 and m 2 , indirect information on its behavior in these 
regions, in the chiral limit, is provided by the two Weinberg sum rules 

rco poo 

Wl= ds Ap{s) = F<° )2 , W2= ds s Apis) = , (39) 
Jo Jo 
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where is the pion decay constant in the chiral limit, and the sum rule for the EM pion 
mass splitting M 



W3 



ee r ds s In -1 Ap(s) = -Ff^Amf 2 , (40) 
Jo J\ oa 



where Amf' 2 is the pion squared-mass splitting in the chiral limit, Amf' 2 ~ m 2 ± —rri^ , and 
we will take A = 1 GeV, for defmiteness, in what follows (the LHS of Eq. (|40| ) is independent 
of A). An additional constraint is provided by the asymptotic OPE form for Ap(s), 

Ap(s) ~ i [BMOi), + B 8 (a 2 s 8 ),] +... , (41) 

valid for sufficiently large s, say s > sa- 

Because of the asymptotic constraint of Eq. (|4"ID , uncertainties in the evaluation of any 
spectral integral are dominated by contributions from the region s ~ 2 — 2.5 GeV — > sa, for 
which spectral data is absent, or has large errors. Letting K(s,p) stand for either Ki(s,p) 
or K 8 (s, p), an obvious way to take advantage of the constraints provided by the chiral sum 
rules is to write K(s, p) in the form 

K(s,p)=C(s,p) + AK(s,p) , (42) 

where C(s,p) is an arbitrary linear combination of the weights occurring in the chiral sum 
rules Wl, W2 and W3, 

C(s,p) =x + ys + zs In (j^ J , (43) 

and the 'residual weight', AK(s,p), is defined trivially by 

AK(s, p) = K(s, p) - C(s, p) . (44) 

The residual weight representation for K(s,p) generates analogous representations for the 
integrals I = 

= ichiral(^) + AJ(^) , (45) 

where 

roc 

/chiral(/i) = / ds C(S, p)Ap(s) = 

Jo 



x-z(-) Ami ) 2 
3a 



(46) 



and 



AI(p)= dsAK(s,p)Ap(s) . (47) 
Jo 

Errors on the determination of /^(/-O associated with uncertainties in our knowledge of 
Ap(s) can then be reduced by adjusting the free parameters x, y, z so as to make AK(s, p) 
small in the region between ~ 2.5 GeV 2 and sa- We will refer to this region as the "matching 
region" in what follows. 
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In order to choose x,y,z in such a way as to minimize the total errors on i^g, it is 
necessary to take into account the uncertainties in our knowledge of the chiral limit values, 
and Amf 2 , which enter the Wl and W3 chiral sum rules. Taking as input the values [17 



Fl 0) = (0.0871 ± 0.0026) GeV and Amf )2 = (0.001174 ± 0.000055) GeV 2 , (48) 

the error, Chiral, on i c hirai(p) becomes 

£ chira i = ±0.000453 GeV 2 \x - 0.674 GeV 2 z\ ± 0.000240 GeV 4 z . (49) 

Based on the NNLO chiral expansion for F n given in Ref. [|17]], we believe that the difference 
between the physical and NNLO chiral values represents an extremely conservative estimate 
of the uncertainty on F(°\ We have therefore taken half this value (±0.0026) as the error 
cited above in an attempt to employ a 'one-sigma' theoretical error akin to that used in 
assigning experimental error. 

The corresponding x, y, z-dependent error on Ai(/i) is obtained as follows. We first 
partition the range of integration into three intervals, 

Ai = [A/] data + [A/] int + [A/| .ay > 

ds AK{s,fi) Ap(s) , (50) 



SA 

+ 



where, as above, sa represents the point beyond which one can employ the OPE represen- 
tation of Ap(s). Our results are insensitive to the actual choice of sa but for definiteness we 
work with sa = 5 GeV 2 . The three regions ('data', 'intermediate' and 'asymptotic') are iden- 
tified by the different ways in which the spectral function Ap(s) is treated. For [Ai]d a ta, the 
compilation of Ap(s) provided by ALEPH is used. In practice, we sum over those experimen- 
tal bins covering the range s < 3.15 GeV 2 . For [Ai] int , we assume |Ap(s)| < Ap max = 0.005 
throughout the interval. Since the three known peaks in Ap(s) decrease in magnitude with 
increasing s (as expected given the 1/s 3 asymptotic fall-off) and since Ap(s) ~ 0.005 for 
s ~ m 2 , this seems to us a reasonable bound. The treatment of [Ai]; nt is discussed in more 
detail below. For [Ai] asy , we take the form given in Eq. fl4T| ) but with a numerical value as- 
sumed for the numerator, viz Ap(s) ~ 7.3 • 10~ 5 GeV 6 /s 3 , which is compatible with our final 
results. The error on AI(p) is obtained by adding the errors associated with each of these 
three contributions in quadrature. The determination of the individual error contributions 
is discussed below. 

The source of the uncertainty on [Ai]<j a ta is Ap(s), the most important component being 
that due to the experimental errors determined by ALEPH as part of their extraction of 
Ap(s) from hadronic r decay data. For each weighted spectral integral, this error is computed 
using the ALEPH covariance matrix. A much less significant uncertainty arises from the 
fact that our theoretical analysis corresponds to the chiral world, whereas our data sample 
is taken in the physical world. It is argued in Ref. [|1|] that such an effect is anticipated to 
be small, (9(m 2 /m 2 ). This expectation is borne out in the analysis of Ref. |19| . 

For the intermediate contribution [Ai] int , we adopt the strategy of assigning the value 

[Ai] int = ± E int , (51) 
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where E- mt is determined by using the Cauchy-Schwarz inequality, 



\AI- mt \ < E- mt = \A Pmax \Js A -m 2 T - J f SA ds [AK{s,fi)} 2 . (52) 

We recall that the quite reasonable bound |Ap max | = 0.005 is employed throughout the 
intermediate region. 

For the asymptotic contribution [AJ] asy , we assume a 100% uncertainty. Although this 
assignment of error is arbitrary, the asymptotic contribution turns out to be so tiny that 
the 100% error does not affect the overall error-minimization in any way. 

At this point, we have a residual weight representation for the integral I (I\ or /§) 
parameterized in terms of three constants x, y, z. Since -E c hirai is independent of y, we 
may reduce E int without adversely affecting -E c hi ra h f° r an Y x , z -, by choosing y so as to 
appropriately minimize \AK(s, fi)\ over the matching region. The lower edge of this region 
is chosen to be so = 2.5 GeV 2 , rather than m 2 ., in order that contributions from that 
portion of the spectrum where experimental errors are large will be suppressed once one has 
performed this minimization. We thus fix y = y(x, z) by minimizing 

(53) 

Like -Echirai, the errors on the three contributions to AI({i) now depend only on x and z, so 
we may combine all errors in quadrature, and minimize the total error with respect to x, z. 
The results of the optimized versions of the ii,s(/-0 residual weight analyses are 



fj, GeV 


h(fi) (GeV 6 ) 


J 8 (/i) (GeV 6 ) 


2 


-(39.7 ±3.1) ■ 


10" 4 


-(26.2 ±3.0) • 


10" 4 


3 


-(63.1 ±5.9) ■ 


10~ 4 


-(31.4 ±5.9) • 


10" 4 


4 


-(82.1 ±9.3) ■ 


lO" 4 


-(34.4 ±9.4) ■ 


10" 4 



(54) 



The uncertainty in each case is acceptable. 

Let us study the various contributions to the integrals at some scale, say, \i = 4 GeV. 
For each integral J, there will be four contributing sources: J c hirai, AJ data , A/ interm and 
A/ asymp . Each source will have an associated uncertainty, except for the chiral contribution 
where we list two errors, the first from and the second from Am[ ' 2 . The numerics are 
as follows: 



Source 


h (4) (GeV 6 ) 


J 8 (4) (GeV 6 ) 


Chiral 


-(23.9 ± 1.4 ±9.1) • 10" 4 


(24.3 ± 1.5 ± 7.5) -10" 4 


Data 


-(58.1 ±0.9) • lO" 4 


-(58.8 ±3.2) • 10~ 4 


Intermed. 


(0±0.7) • 10~ 4 


(0±4.6) • 10" 4 


Asympt. 


-(0.1 ±0.1) • 10~ 4 


(0.1 ±0.1) • io- 4 



The largest source of uncertainty in our results turns out to originate with the 4.6% er- 
ror in Eq. (|48| ) for the pion squared-mass splitting (equivalent to Amp D ' = 0.32 + / — 
0.20 MeV). Our knowledge of this (m u - m d ) 2 effect arises mostly from studying the 
kaon mass difference. There, the effect of electromagnetism (which shows significant SU(3) 



ds 



-S[| 



K (s, fi) — x — ys — zs In 



A 2 
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breaking) is subtracted and the remaining kaon mass difference fit by — m u in lowest 
order chiral perturbation theory. However, the next order chiral corrections could easily be 
sizeable at the kaon mass scale, and one must also cope with ir -(rf,rf') mixing. In light of 



these effects, we feel the current 4.6% error cited in Ref. |l7j is realistic. 

Finally, we point out that a possible modification of the procedure described above 
would be to include the F^ 3 factors of Eq. (|3|) ab initio in our error minimization analysis. 
In other words, we would be determining integrals J^g = h^/F^ 3 rather than just Z^g. 
Such a program produces results completely compatible with those of Eq. ([54]) and with no 
significant reduction in the size of the errors. 



B. Renormalization Group Evolution and Results 

In this section we outline the procedure for arriving at the physical matrix elements 
((7T7r)/ =2 | Q 78 |i^ ) in the chiral limit. 

Using the values for Ii,g(/i = 4 GeV) obtained in the previous section and setting Hi t s(fi = 
4 GeV) = 0, we solve Eq. ( |3"3"D for (01,8)^=4 GeV- Since we work at the scale \i = 4 GeV we 
use the values of our coefficients corresponding to four active flavors (jif = 4). Moreover 
we use ^qcd = (340 ± 50) MeV. We subsequently apply the NLO renormalization group 
evolution, to obtain (01^)^=2 GeV- Denoting by U(^l,^h) the operator governing the NLO 
evolution of the Wilson coefficients (ci(/z), c 8 (/z)) = ^(aO between the scales /jl and fin, 

c(/i L ) = U(fi L ,fJL H ) c(fi H ) , (56) 

the local operators (Oi(fi), 8 (//)) = T (fi) evolve according to: 



(O)M ■ (57) 



The results for (C^g)^ and ((7r7r)/ =2 | Qi^K®)^ are reported in Table [TJ for the scales \i = 
2, 4 GeV. The vacuum matrix elements (Oi, 8 )^ are converted to the K — > tttc matrix elements 
through Eqs. ([3D. We have added in quadrature the uncertainties from the vacuum matrix 
elements and from the 1 / F^ 3 factor. 



IV. FINAL COMMENTS 

Our goals in this paper have been to exactly match the dispersive description of the 
weak amplitudes {(7T7r)7 = 2| Qt.sI-^" ) i n the chiral limit to the specific conventions of the 
existing two-loop calculations of the Wilson coefficients which appear in the effective weak 
hamiltonian, and to provide a numerical evaluation using as input only experimental data 
and the rigorous theoretical constraints of QCD. 

In order to accomplish the matching, we adopted the operator basis and renormalization 
conventions of Ref. 0. On one hand we analyze the OPE by matching the hamiltonian at 
the scale Mw and evolving the operators to a lower scale fx. By simultaneously considering 
the dispersive representation of the matrix element, we use the OPE to identify the MS 



13 



TABLE II. Numerical values at /x = 2 GeV and \x = 4 GeV. 





(C 1 ) M (GeV 6 ) 




(GeV 6 ) 


H (GeV) 


NDR HV 


NDR 


HV 


2 


-(0.53 ± 0.34) • 10" 4 -(1.64 ± 0.18) • 10" 4 


-(14.4 ±4.3) • 10" 4 


-(15.2 ±4.4) • 10" 4 


4 


-(1.47 ± 0.18) • 10" 4 -(2.59 ± 0.33) • 10" 4 


-(19.3 ±5.4) • 10" 4 


-(20.3 ±5.6) • 10~ 4 




((vrvr) /=2 |Q 7 |^°) M (GeV 3 ) 


((vr7r)/ =2 |Q8 


(GeV 3 ) 


H (GeV) 


NDR HV 


NDR 


HV 


2 


0.16 ±0.10 0.49 ±0.07 


2.22 ±0.67 


2.46 ± 0.70 


4 


0.44 ±0.07 0.78 ±0.12 


3.06 ± 0.87 


3.32 ± 0.90 



operators at the scale \i in both the NDR and HV schemes. Our results also identify the role 
of higher dimensional operators, which enter the dimensionally regularized matrix elements. 

In the numerical analysis we have developed the 'residual weight' method, which provides 
an evaluation requiring as input only data and the chiral-limit values and Amjf'. This 
method expresses the dispersive integrals Ii(fi) and Is(fi) in terms of known quantities (the 
other known dispersive sum rules) plus a residual integral. For the residual integral, we make 
use of the ALEPH data and employ only a conservative bound on the spectral function in 
the region where r decay data is not available. We have chosen to evaluate Ii(fi) and /s(/-0 
at large values of \x where the effects of higher dimension operators are negligible. We then 
use the renormalization group to obtain the values of the matrix elements at other scales 

Our results lead to the matrix elements given in Table II and Table III. Because we only 
use real experimental input, in combination with very conservative assumptions about the 
uncertainties on Fj°) and Amf, the final errors we quote are very conservative ones. These 
are the maximal error bars of the dispersive evaluation — our method has been designed to 
give as pure and conservative an evaluation from data as possible, and other methods may 
reduce the error bars. For example, Eq. (Q) shows that the uncertainty on the residual 
weight analysis decreases dramatically at lower values of //. To directly work at // = 2 GeV 
one needs to evaluate, and correct for, the effects of higher dimension operators. We have 
been studying this problem through the use of FESR and preliminary indications yield 



results consistent with the present evaluation but with reduced error bars [T5|]. However 
that analysis requires the development of further techniques and consistency checks, and 
we will present it in detail elsewhere ||16|| . The present method stands on its own as being 
especially simple and model independent. 

In Table |T|, we compare our results at the scale fj, = 2 GeV to other determinations: 
the lattice Monte Carlo simulation [2~u], the large N c approach |I2|], the so-called X-boson 
approach |2]J and the vacuum saturation approximation (VSA). We make the following 
remarks regarding this comparison. 

The Monte Carlo simulation of Ref. pOl is displayed in the 'Lattice' row of Table III]. 
Note that their HV results have been converted to the HV scheme used in this work. As 
regards comparison between the lattice results and ours, we should point out that the two 
approaches correspond to somewhat different limits. The lattice evaluation extrapolates the 
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results to the physical values of the quark masses, while our dispersive analysis applies to 
the amplitude in the chiral limit. Even so, the extrapolation between the chiral limit and 
the physical masses has been addressed by both analytic p2lp3j and lattice methods |24 



and the differences are not large enough to account for the discrepancy. This means that the 
difference between the lattice and dispersive values for Q$ represents a serious disagreement. 
It will be important to pursue an understanding of the physics behind this disagreement, as 
in principle both methods are fully rigorous. 

In the case of the large N c evaluation [|1^], there are differences in principle with our 
approach because they have not performed the complete matching at two loops. In addi- 
tion, instead of using data for the evaluation of the vacuum polarization functions, they 
use a vector meson dominance model. Nevertheless, their results are consistent with ours, 
especially once the NLO radiative corrections are included in their work. 

A recent manuscript also discusses matching and dimension-eight effects, but in the 



context of an 'X-boson method'. The greatest procedural difference with our work comes in 
the evaluation of the Q 8 matrix element. They find small non-factorizeable effects such that 
their dominent contributions are the factorized terms, refined to include the scale dependence 
from the matching procedure. 

The VSA values are obtained from 

IF 

((^)i=2\Qs\K°}y SA = S((7ttt) i=2 \Q 7 \K°} vsa = — ~ 0.94 GeV 3 , (58) 

[m a + m d )^ =2 GeV 

where we take {m s + md)^=2 GeV = HO MeV based on the value of m s given by the midpoint 
of the combined sum rule and lattice ranges in the recent review Ref. |25| in combination 
with quark mass ratios obtained in ChPT analyses 



The matrix element for Q 7 has little phenomenological interest, and will be useful pri- 
marily for the comparison with the lattice as both the dispersive and lattice evaluations 
improve. However the matrix element for Qg is one of the main contributions to e'/e, and 
gives a negative contribution to that quantity. Our rather large result implies that the 
electroweak penguin contribution is 



= (-2.2 ±0.7) x 10~ 3 . (59) 

EWP 



The gluonic penguin matrix element needs to be large enough to bring the result up to the 
experimental value. Given this important consequence, the resolution of the disagreement 
of the lattice and dispersive values will be specially interesting. 
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TABLE III. Comparison of Matrix Element Determinations at /i = 2 GeV. 



<(^r) /=2 |g 7 |^°>(GeV 3 ) ((^vr) /=2 |Q 8 |E: )(GeV 3 ) 



Method of Calculation 


NDR 


HV 


NDR 


HV 


This work 
Lattice @ 
Large N c |12| 
X-boson |21| 
VSA 


0.16 ±0.10 
0.11 ±0.04 
0.11 ±0.03 
0.26 ±0.03 
0.32 


0.49 ± 0.07 
0.18 ±0.06 
0.67 ±0.20 
0.39 ±0.06 
0.32 


2.22 ±0.67 
0.51 ±0.10 
3.5 ± 1.1 
1.2 ±0.5 
0.94 


2.46 ± 0.70 
0.57 ±0.12 
3.5 ± 1.1 
1.3 ±0.6 
0.94 
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